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Abstract 

The mean-field theory represents a useful starting point for studying carrier-induced ferromagnetism in Mna;IIIi_a; V 
diluted magnetic semiconductors. A detail description of these systems requires to include correlations in the 
many-body hole system. We discuss the effects of correlations among itinerant carriers on magnetic properties of 
bulk Mna;IIIi_2,V and magnetic semiconductor quantum wells. Presented results were obtained using parabolic 
band approximation and we also derive a many-body perturbation technique that allows to account for hole-hole 
correlations in realistic semiconductor valence bands. 
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1. Introduction 

Recent advances [1-3] in the fabrication and con- 
trol of Mna;IIIi_a;V diluted magnetic semiconduc- 
tors (DMS) have opened up a broad and relatively 
unexplored frontier for both basic and applied re- 
search. We have developed a mean-field theory of 
free-carrier induced ferromagnetism in DMS which 
is intended to be useful for bulk materials as well 
as for any spatially inhomogeneous systems. Spe- 
cial emphasis is placed on interaction efl^ects in the 
itinerant many-body system whose inclusion is re- 
quired for a detailed understanding of magnetic 
properties of these materials. In Section 2 we derive 
the mean-field theory equations. Correlation ef- 
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fects on magnetic properties of bulk DMS and Mn- 
doped semiconductor quantum wells are studied 
in Sections 3 and 4, respectively, using a parabolic 
band model. In Section 5 we present a many-body 
perturbation theory that accounts for hole- hole in- 
teractions in realistic semiconductor valence bands 
and evaluate the exchange enhancement of the fer- 
romagnetic critical temperature for typical sample 
parameters. 



2. Mean-field theory 

Our theory is based on an envelope function de- 
scription of the valence band electrons and a spin 
representation of their kinetic-exchange interac- 
tion [4] with d-electrons on the S = 5/2 Mn++ 
ions: 
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H = Hrr,+Hf + JpdJ2^I- - ^l)^ (1) 

i,I 

where i labels a free carrier, / labels a magnetic 
ion and the exchange interaction is parametrized 
by a coupling constant Jpd- In Eq. (1) Hm is the 
Hamiltonian of the magnetic ions. Hj is the six- 
band Luttinger Hamiltonian for free carriers in the 
valence band, S is the magnetic ion spin and s is 
the electron-spin operator projected onto the j = 
3/2 and 1/2 valence band manifold of the Luttinger 
Hamiltonian. 

In the absence of external fields, the mean po- 
larization of a magnetic ion is given by [5] 

{m)i = SBs{JpdS[n^{Iii) - n^CR/)] /2fcBT) , (2) 

where Bs{x) is the Brillouin function, 

Bs{x)^^^x, x^l. (3) 
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The electron spin-densities no-(r) are determined 
by solving the Schrodinger equation for electrons 
which experience an electrostatic potential, Wes(r), 
a spin-dependent kinetic-exchange potential. 



hpd{r) = Jpd^(5(r-R7)(m)j 



(4) 



and a local-spin-density-approximation (LSDA) 
exchange-correlation potential, Wa;c(r), on which 
we comment further below. The kinetic-exchange 
potential is non-zero only in the ferromagnetic 
state and we assume that the magnetic ions are 
randomly distributed and dense on a scale set by 
the free carrier Fermi wavevector and that their 
density, c(r), can be controlled. This allows us to 
take a continuum limit where 



400 



300 




n (10 cm ) 



Fig. 1. Curie- Weiss temperature as a function of hole den- 
sity calculated including (solid line) and neglecting (dashed 
line) hole exchange-correlation potential. Inset: relative 
spin polarization of holes as a function of temperature for 

n = 10^" cm~^. 

where %/ is the interacting hole magnetic suscep- 
tibility. The Curie- Weiss temperature, obtained 
from Eqs. (2) - (6), is given by 



kBTc = 



cSjS + 1) 4dXf 
3 (5MB )^ 



(7) 



hpd{r) = Jpd c(r) (m)(r). 



(5) 



In Fig. 1 we have plotted the ferromagnetic tran- 
sition temperature as a function of the total car- 
rier density, n, predicted by this expression for p- 

type Mnj;Gai_2;As with Jpd = 0.15 cV nm^ c = 
lO^i cm"*^ and for valence bands approximated by 
a single parabolic band with effective hole mass 
m* = CSrUe. When exchange and correlation ef- 
fects are neglected and x/ is replaced by its zero- 
temperature value, Tc is proportional to n^^^. In- 
cluding the parabolic band exchange-correlation 
potential [6] enhances Tc by » 20% for typical hole 
densities 



3. Homogeneous DMS 



4. DMS quantum wells 



For homogeneous systems with randomly dis- 
tributed localized spins the mean-field equations 
can be solved analytically. The hole spin-density 
and the kinetic-exchange potential are related as 

Xf 



(n| - ni)/2 = 



'pd 



(6) 



In the parabolic band approximation, the spin- 
densities in Mnj;Gai_2;As quantum wells arc de- 
termined by solving the Schrodinger equation 



2m| ^ dx'^ dy"^ ' 



1 92 
2mt dz^ 



2 



+Vxc,a{z) - -^hpd{z) 



(8) 



where the in-plane effective mass m| « 0.1 Imo 
and the out-of-plane mass m* w 0.38mo. The elec- 
trostatic potential varies along the growth {z) di- 
rection and includes band offset and ionized im- 
purity contributions. The LSDA equation for the 
exchange-correlation potential reads 



d[nexc{n^, ni)] 
drier 



(9) 



where e^^cl^^T'^^i) exchange and correlation 

energy per particle of a spatially uniform system 
[6]. 

If hole correlations and subband mixing are ne- 
glected the following anal}d;ic expression can be 
found for Ty. 



S{S + 1) 4d ^\\ 
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dz\ip{z)fc{z) . (10) 



In Eq. (10), the critical temperature increases 
with the subband wavefunction, ip{z), and mag- 
netic impurities overlap, is proportional to J^^ 
and m|. Similar Tg-equation was derived [7] us- 
ing the RKKY theory and its predictions are in 
agreement with the measured carrier-induced fer- 
romagnetism in II- VI quantum wells [7,8] . Here we 
emphasize the role of hole-hole interactions and 
Mn doping profile on magnetic properties of biased 
III-V DMS quantum wells. A remarkable feature 
of quasi-2D ferromagnetic DMS is the possibility 
of tuning Tc through a wide range in situ, by the 
application of a gate voltage. In Fig. 2 we illustrate 
this for the case of a w = 10 nm quantum well with 
magnetic ions covering only a d = 3 nm portion 
near one edge. Even for such a narrow quantum 
well, the critical temperature can be varied over 
an order of magnitude by applying a bias voltage 
which draws electrons into the magnetic ion re- 
gion. Hole-hole interactions substantially increase 
the critical temperature, as also shown in Fig. 2 
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Fig. 2. Dependence of Tc on bias voltage applied across 
the w = V) nm wide quantum well partially occupied by 
magnetic ions over a distance of 3 nm near one edge of the 
quantum well. Circles (squares) represent results of the full 
numerical self-consistent calculations without (with) the 
exchange-correlation potential. 

5. Excheinge cind correlations in realistic 
valence bands 



A semiempirical six band model is used for quan- 
titative description of electronic states in cubic 
semiconductors near the top of the valence band. 
For systems with spin-orbit coupling the Kohn- 
Luttinger Hamiltonian [9] is derived in the basis of 
total angular momentum eigenstates with j = 3/2 
and 1/2, and = — j, —j + 1, j. The many- 
body interaction Hamiltonian in this basis can be 
written as 

^1= "i2,k2-q«Il,ki-Kqani,kian2,k2V'(q),(ll) 

ni,''i2, 
ki,k2,q 

where V{q) is the Fourier transform of the 
Coulomb interaction. Using the following notation 
for the transformation between creation and an- 
nihilation operators for the total angular momen- 
tum eigenstates and Kohn-Luttinger Hamiltonian 
eigenstates: 

^n,)s. ~ ^ ] Zn^a^k ^a.ki <ln,k = ^n,a,k Q^a,k (12) 

a a. 

the exchange contribution to the hole total energy 
has a form 



3 



ai,a2, 
k,q 



E 



'2'n,a2,k+q^n,ai,k 



nq). (13) 



where na,k is the Fermi distribution function of 
Kohn-Luttinger eigenstates. 

In the random phase approximation (RPA) , the 
correlation energy is related to the wavevector and 
frequency dependent dielectric function as [10] 



cko 
2^ 



tan 



lmeRPA{q,u)) 



R.eeKPA(q,w) 



eRPA{q,tu) = l-y(q)i3°(q,a;) , 



(14) 



where -D"(q, uj) is the semiconductor valence band 
density-density response function: 

-D°(q,a;)= ^ nai,k(l - rias.k+q) 

ai ,Q2, 
k 



^ ^ -^71,02 ik+q-^^Tijai ,k 

1 



+ (Cai,k - ^c«2,k+q) + iV 
1 



?i<^ - (^ai,k - Ca2,k+q) - iV . 



(15) 



and ^a,k arc the Kohn-Luttinger eigenenergies. 
The Curie- Weiss temperature for realistic valence 
bands is obtained from Eq. (7) with the magnetic 
susceptibility given by 



Xf 



(16) 



where the itinerant system total energy per vol- 
ume, Et/^, has a kinetic (band) contribution and 
the exchange and correlation contributions given 
by Eqs. (13) and (14). 

For the current experimental values of the hole 
density n 
c = 1 X 10- 
stant Jpd = 55 meV nm^, the six-band model 
gives mean-field critical temperature for bulk 



= 3.5 X 10^° cm ^, Mn concentration 
"•^^ cm~^, and kinetic exchange con- 



Mna;Gai_xAs, Tc = 86 K when exchange and 
correlations are neglected. The exchange and cor- 
relation contributions to the total energy enhances 
Tc by a factor of 1.2. Compared to the experimen- 
tal Tj. = 110 K, the mean-field theory accurately 
predicts the ferromagnetic critical temperature of 
Mna;Gai_a;As DMS. 
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